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Abstract 

We study the sensitivity of presupernova evolution and supernova nucleosyn- 
thesis yields of massive stars to variations of the helium-burning reaction rates 
within the range of their uncertainties. The current solar abundances from Lod- 
ders (2009) are used for the initial stellar composition. We compute a grid of 
12 initial stellar masses and 176 models per stellar mass to explore the effects 
of independently varying the 12 C(o;,7) 16 and 3a reaction rates, denoted R a ,i2 
and i? 3ct , respectively. The production factors of both the intermediate-mass el- 
ements (A=16-40) and the s-only isotopes along the weak s-process path ( 70 Ge, 
76 Se, 80 Kr, 82 Kr, 86 Sr, and 87 Sr) were found to be in reasonable agreement with 
predictions for variations of R 3a and R a ,i2 of ±25 %; the s-only isotopes, how- 
ever, tend to favor higher values of R 3a than the intermediate-mass isotopes. 
The experimental uncertainty (one standard deviation) in Rz a {R a ,i2) is approx- 
imately ±10%(±25%). The results show that a more accurate measurement of 
one of these rates would decrease the uncertainty in the other as inferred from 
the present calculations. We also observe sharp changes in production factors 
and standard deviations for small changes in the reaction rates, due to differ- 
ences in the convection structure of the star. The compactness parameter was 
used to assess which models would likely explode as successful supernovae, and 
hence contribute explosive nucleosynthesis yields. We also provide the approxi- 
mate remnant masses for each model and the carbon mass fractions at the end 
of core-helium burning as a key parameter for later evolution stages. 

Subject headings: Nuclear Reactions, Nucleosynthesis, Abundances, Sun: Abundances, 
Stars: Supernovae: General 
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Introduction 



Massive stars are responsible for the production of most intermediate- mass (A = 16—40 ) 



isotopes through hydrosta tic burning phases and subsequent supernovae (IBurbidge et al 



1957 



Woosley et al 



2002f ). During core-He burning the 12 C(a,7) 16 and 3a reactions 
compete to determine the relative abundances of oxygen and carbon prior to core-C burning. 
Changes in these abundances have significant effects on subsequent stellar evolution and 
structure and on the resulting nucleosynthesis. The carbon abundance influences subsequent 
shell burning episodes and affects whether core-C burning will be radiative or convective. 
There is also a non-monotonic rela tion between the carbon abundance and the resulting 



remnant mass (jWoosley et al 



20031 ). so that these reactions are important for understanding 



the populations of neutron stars (NS) and black holes (BH). 

These rates can also affect weak s-process yields. The weak s-process is a slow neutron 
capture process occurring at th e end of convective core-He burning and during shell carbon 



burning (jPignatari et al. 



2010), and contributes to isot opic production alo ng the s-process 



19931 ). A change in 



path above iron and up to a mass number of A^lOO (IRaiteri et al. 
the helium-burning rates can induce a corresponding change in temperature to keep the 
star at constant luminosity, and the reaction for the neutron source for the weak s-process, 



22 T 



Ne(a,n) Mg, is highly temperature-dependent. Additionally, the amounts of neutron 



poiso ns have been shown to vary with these rates ( iRayet and Hashimoto 
200% . 



2000 



Tur et al 



Although not dis cussed in 
26 Al, 44 Ti, and 60 Fe (ITur et al 



;his p aper the production of the important radioactive nuclei 



20101 ) is also sensitive to these rates. Gamma rays from 
these nuclei provide observational information that may help to test models of massive star 
internal structure and nucleosynth esis through the c onstraints imposed by t he abundance 



ratios of 44 Ti/ 56 Co and 26 Al/ 60 Fe flDiehl et all 



2006 



Leising and Diehl 



20091 ) 
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This work is an extension of the study by Tur et al. (120071 . 120091 ). who calculated a 
limited subset of models of our full 2D parameter space. They concluded that across the 2a 
uncertainty range root-mean-square (rms) deviations for the production factors sometimes 
vary non-monotonically with the rates, as do deviations in the re mnant mass, ind icating; 



2003). We 



that both helium burning reactions are independently important (jTur et al. 
extend upon their study by performing a much finer sampling of the 2a uncertainty range 
to assess the effect of changing the rates independently, in order to map the whole 2D 
landscape to determine the true non-monotonic behavior with a higher resolution grid. The 
purpose is to: i) examine the effect of varying the helium burning reaction rates on the 
production factors of intermediate-mass isotopes, ii) examine the effect on the production 
factors of the six s-only isotopes along the weak s-process path, iii) assess the impact on i) 
and ii) of including only models that are likely to explode as successful supernovae, and 
iv) examine the effect on the remnant mass. We do not address the effect of varying the 
solar abunda nce set, w h ich in deed has been shown to impact the final n ucleosynthesis, an d 



is studied in 



Tur et al. 



Tur et al. 



2009) 



2009|) 



(120071 ) for the intermediate-mass isotopes and in 
for the weak s-process isotopes. We use an updated solar abundance set ( lLodders 
that was not yet available for the previous studies. It was corrected for the weak s-only 
isotopes by subtracting estimated main s-process contributions, as described in Section 2. 



We note that recent 3-body calcul ations of the c 
at temperatures below ~ 0.07 GK (jNguyen et al. 



a reaction show an increase in this rate 



20121 ). This will not impact the current 



study as He ignition in massive stars occurs beyond this threshold, at T 3q ~ 0.1 GK. 

This paper has the following outline: in Section 2 we describe the stellar models 
and methods used in the analysis. In Section 3 we compare the intermediate-mass and 
weak s-process isotopes across all models, and for the subset of models likely to explode 
as supernovae rather than collapse to black holes (we ignore hypernovae and gamma ray 
bursts). This subset is chosen using a compactness parameter filter. We also discuss the 



-5 - 



remnant mass and carbon mass fractions at the end of core-He burning for the different 
stellar masses. Our conclusions are given in Section 4. 



2. Stellar Models and Analysis 



All models were computed using KEPL 



hydrodynamics package for stellar evolution (IWeaver et al- 



ii, a time-implici t one-dimensiona l 



1978 



Rauscher et al. 



20021). 



A grid of 12 initial stellar masses M/M =12, 13, 14, 15, 16, 17, 18, 20 22, 25, 27, and 



30 was used, with the revised pre-solar abundances from 



Lodderd (120091 ) for the initial 



composition. For each stellar mass, 176 models were computed to scan at least a 2a 
uncertainty range for both the 12 C(a,7) 16 and 3a reaction rates, denoted R a ^2 and R, 



3«; 



respectively. This 



rates, as done by 



range wa s para metrized as a multiplier on the centroid values for the 



Tur et al. 



(120071 ). The centroid values used in KEPLER are taken from 



Caughlan and Fowler! ([1988) for i? 3a , and 1.2 times the rate recommended by 



Buchmann 



( 20001 ) for R a) \2- The range for the R a ^2 multipliers was (0.5,2.0) with a resolution of 



A = 0.1, and the range for the R^ a multipliers was (0.75, 1.25) wi th a resolution o 



Commonly accepted uncertainties for R% a and R Q) \2 are ±10 % (jChernykh et al. 



A = 0.05. 



20101) and 



±25 %. Our total range was, conservatively, slightly more than ±2er. 



In Fig.[T]we show the grid of models computed both by 
present work. 



Tur et al. 



(120071 ) and in the 
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Fig. 1. — The R at u and i?3„ reaction rate multiplier values used in the stellar models. The 
models using the reaction rate multip l ier pa irs given in red plusses, green asterisks, and blue 



Tur et al. 



(120071 ) . The black asterisks show the models computed 



squares were performed by 
in the present work. 

All stellar models were first evolved through hydrostatic burning until the Fe core 
collapsed, and an inward velocity of 10 8 cms _1 was reached. The explosion mechanism for 
the resulting supernova was modeled as a mechanical piston that imparted an acceleration 
at constant Lagrangian mass coordinate to provide the desired total kinetic energy of the 
ejecta, taken in these models to be 1.2 B (IB = 10 51 erg) at 1 year aft er the explosion. Fo r 



detai 



see 



s on the parametrization of the explosion used in KEPLER see IWoosley and Heger 



(120071). and references therein. For details on the treatm ent of convection and mixing 



Woosley and Weaver 



cut is given in 



Tur et al 



(1998 



) and 



(120070 and 



Wooslev et al. 


( 


2002 


), and 


Heeer and Wooslev 


(2010) 



, and a discussion of the mass 



magnetic fields are ignored. The final supernova yields for all models were then averaged by 
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integrating over the Salpeter initial mass function (IMF) for each reaction rate multiplier 
pair in Fig.[TJ The yields from stellar winds were included. 



In Equation[T] the IMF interpolated yield mass for isotope i is given by Y*, and the 
Salpeter mass spectrum is £ (m) = C -m -2 ' 35 , where C is the proportionality constant . The 
mass grid used for the integrations are the ejected masses, defined as the baryonic remnant 
masses subtracted from the initial stellar mass grid. Yields are linearly interpolated between 
adjacent masses in the integral, with the slope defined by, Sij = — Y it j) / ( 

where Yij is the yield mass of isotope i from a model with initial mass rrij. In Equation[2] 
the production factor for isotope i is given by Pi, the sum in the denominator runs over all 
isotopes, and the solar mass fraction of isotope i is given by Xj i0 . 

Massive stars are responsible for producing nearly the entire solar abundance of a subset 
of intermediate-mass isotopes, namely 16 ' 18 0, 20 Ne, 23 Na, 24 Mg, 27 Al, 28 Si, 32 S, 36 Ar, and 
40 Ca. Hence, in order to make the solar abundance pattern, massive near-solar metallicity 
stars need to produce these intermediate-mass isotopes in solar ratios. An analysis of the 
standard deviations of the production factors for this set of isotopes is used in the present 
work to identify helium reaction rate values that agree with solar observations. This 
agreement is an approximation that the above isotopes owe their entire solar abundance 
to massive, near-solar metallicity stars, and relies on sufficient sampling of the initial mass 
function (IMF) and understanding of the initial compositions. As mentioned, the impact of 
uncertainties in the initial composition is not addressed in this work, but an analysis of the 






(2) 
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effect of different compositions can be found in iTur et al.l ( 120071 ) 



For the weak s-process isotopes a correction is necessary, since the solar abundances 
of the six s-only isotopes along the weak s-process path ( 70 Ge, 76 Se, 80 Kr, 82 Kr, 86 Sr, and 
87 Sr) have additional contributions from the main s-process, which occurs in asymptotic 
giant branch (AGB) stars, not massive stars. Hence what is needed are the production 
factors relative to the contributions from the weak s-process only, not relat ive to the entire 



solar a bundance. To achieve this, the solar abundance decomposition from 



West and Heger 



( 120121 ) was used, which gives, in part, the approximate solar contributions for the six weak 
s-only isotopes. This modifies Equation[TJ with the substitution of X i<w for X ijQ , where X i<w 
denotes the contribution to the solar mass fraction of isotope i from the weak s-process. 



3. Results and Discussion 

3.1. Comparison of 12 C(a,7) ls O and 3a Reaction Rates 

The production factors for the models were computed using Equations!]] and 121 and 
standard deviations were calculated for each model, using the intermediate-mass isotope set 
(in Section 4.1) and the six s-only isotopes along the weak s-process path (in Section 4.2), 

where are the production factors computed in Equation[2], and n is the number of entries 
in the isotope list: 10 for the intermediate-mass isotopes and 6 for the s-only isotopes along 
the weak s-process path. We distinguish op (the standard deviation of the production 
factors) from a (the uncertainty in the rates). Since massive stars contribute to most of 
the solar abundances for the intermediate-mass isotopes considered (or a fraction of them 
in the case of the s-only isotopes), low standard deviations should indicate combinations of 
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R n 12 and R 



by 



Tur et al. 



that agree with observations. We first follow the type of analysis performed 
(j2007l ). for constant R a ,u and R 3a multipliers, shown in Fig. [21 
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Fig. 2. — Left: The production factors for the intermediate-mass isotopes averaged over the 
IMF as a function of the -R Qj i2 multiplier, at a constant R^ a multiplier of 1.0. Right: The 
production factors for the intermediate-mass isotopes averaged over the IMF as a function of 
the R 3a multiplier, at a constant R a ,i2 multiplier of 1.2. The inset graph shows the standard 
deviations. 



For a constant R 3a multiplier of 1.0 (Fig. [21 left), values for the R at i2 multiplier are 
favored within about ±25 % of the centroid multiplier of 1.2. For a constant R a ,i2 multiplier 
of 1.2 (Fig. 121 right), values for the R 3a multiplier have a minimum standard deviation at 
0.85, and vary across a similar range of standard deviation values for a ±25 % change in 
R% a - Since R% a is better experimentally determined, however, the extremes of this range are 



less likely than for R, 



the fin dings of 



use the 



aJ2. 



Lodders 



ur et al. 



he re sults in Fig. [2] show approximate qualitative agreement with 



(120071 ). but care must be taken in a comparison as they do not 



20091 ) abundances, and they have demonstrated that there is non-trivial 



variation among different solar abundance sets. 



The corresponding plot for the weak s-only isotopes is given in Fig.O 
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l2 C(a,7) ,G Rote Multiplier (Triple-a constant at 1.0) Triple-a Rote Multiplier ( 12 C(a,7) 1B constant at 1.2) 

Fig. 3. — Left: The production factors for the weak s-only isotopes averaged over the IMF as 
a function of the i? Qj i2 multiplier, at a constant R 3a multiplier of 1.0. Right: The production 
factors for the weak s-only isotopes averaged over the IMF as a function of the R 3a multiplier, 
at a constant R a ,i2 multiplier of 1.2. The inset graph shows the standard deviations. 

For a constant R 3a multiplier of 1.0 (Fig.[3l left), values for the R a ,i2 multiplier have a 
minimum standard deviation at 1.3. For a constant R a ,i2 multiplier of 1.2 (Fig.[3l right), 
the standard deviation has a minimum at the R% a multiplier value of 0.95. Significant 
variations in the production factors exists across both multiplier ranges. 

3.2. Intermediate-Mass Isotopes 

To address how changing the rate multipliers independently affects the nucleosynthesis, 
the entire set of models in (R a ,i2,R3a) space was mapped in a 2D grid, rather than 
restricting ourselves to ID slices. Results are first given for just the 25 M Q models (Fig.Hj). 
The corresponding plots for all models can be found in the Appendix. 
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Fig. 4. — Standard deviations of the production factors as a function of the R a ,x2 and R^ a 
reaction rate multipliers for the 25 M models. Each model is given by the R a ,i2 and R^ a 
reaction rate multiplier pair used for the helium rates. 

The best fit -R Qj i2 and R% a values occupy a region in the lower right-hand corner, and 
a strip running approximately though the centroid value for each rate. It is interesting 
that some adjacent models display significant differences in their nucleosynthesis despite 
small change in the reaction rate multipliers, for example (R a ,i2, R$a) — (0.5, 0.9), (0.6, 0.9) 
or (R a ,i2,R3a) — (0.5, 1.05), (0.6, 1.05). In some cases adjacent models can evolve with 
quite different shell burning episodes, whereas in other cases can be very similar, such as 
for {R a ,n,Rz a ) = (0.9, 0.9), (1.0, 0.9) or (R a ,i2, Rsa) = (1.0, 1.1), (1.1, 1.1). How this occurs 
can be understood by considering the convective history of adjacent models. First, the 
(Ra,i2, R3a) = (0.5,0.9) and (R a ,i2, R3a) = (0.6,0.9) models are given in Fig.[5l which have 
a difference in ap of 3 to 4 (Fig.HJ). 
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6 4 2 0-2-4-6-8 

log( time till core collapse / yr ) 



Fig. 5. — Top: The convection plot for the inner 7M of the (-R a ,i2 ; ^3a) = (0.5, 0.9) model. 
Bottom: The convection plot for the {R a ,i2, Rsa) = (0.6,0.9) model. Shown are convective 
regions (green hatch- marks), semi-convective (red cross-hatching), energy generation from 
nucleosynthesis (blue), and radiative/neutrino cooling (pink). The entire evolution from the 
main sequence to onset of core collapse is shown. 

In the (0.5, 0.9) model we observe a convective region that extends past the C shell and 
into the above He layer, with subsequent He ingestion into the shell burning with C into O. 
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In the (0.6, 0.9) model the convective layer terminates before the He shell, and the mixing 
and subsequent nucleosynthesis seen in the (0.5, 0.9) model does not occur. 

In contrast, the (R a ,i2, Rza) = (0.9,0.9) and (-R a ,i2, R3a) = (1.0,0.9) models are given 
in Fig.0 which show a difference in ap of < 1. 




6 4 2 0-2-4-6-8 

log( time till core collapse / yr ) 



Fig. 6. — Top: The convection plot for the inner 7M Q of the (R a ,i2, R3a) — (0.9, 0.9) model. 
Bottom: The convection plot for the (R a ,i2, Rsa) — (1-0, 0.9) model. See Fig.[5]for a detailed 
description. 
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In both Fig.[5]and Fig. [6] the AR^ a is the same; however, in the latter the R ay u values 
are sufficiently large to result in radiative core-C burning, and C shell burning episodes in 
the (R a ,i2, Rsa) — (0.9,0.9) and (1.0,0.9) models (Fig.EJ) that do not interact convectively 
with the He layer. Generally, larger carbon abundances at the end of core-He burning can 
support longer and more energetic carbon shell burning episodes, which can result in He 
ingestion leading to different nucleosynthesis. 

The differences of adjacent models depend on the initial stellar mass, and it is expected 
that regions in the {R a ,i2iRza) parameter space that have different standard deviations at 
one mass may not at another. To average these effects across the IMF, Equations[TJ and [5] 
were used for the entire set of models for all masses, and the standard deviations for the 
production factors are given Fig. [7J 

Changes in the 12 C + 12 C reaction rate may also affect nucleosynthesis; rate uncertainties 
at low temperatures were though t to be large, perhaps orders of magnitude. A heavy 



ion fusion study 
the other hand, 



by_ 



Jia ng et a 



Spilla ne et al 



(2007) reported a rate decrease at low energy. On 



(2007) found a strong increase in rates due to a low 



energy resonance. A still lower energy resona nce at lower ene r gy ha d weak experimental 



support. Possible effects have been studied by lPignatari et al. 



(120131 ) for a 25 M half-solar 



metallicity star. They found significant changes in the production of s-process elements. 
However, their calculations were not carried thr ough to solar co llapse and an explosion, 



which affects the production of these isotopes 



measurements by 



( ITur et al 



Zickefoosd (12010. unpublished! ) . reported by 



2009[). In addition, re cent 



Notani et al. 



(120121 ) give an 



S-factor 50 times smaller than the previously reported value. For these reasons it is not 
clear how much the nucleosynthesis considered in this paper would be affected. 



- 15 - 
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,2 C(a,7)' 6 Rate Multiplier 



Fig. 7. — Standard deviations of the IMF averaged production factors for the intermediate- 
mass isotopes. The entire grid of initial masses was used. 

The results in Fig. [7] show a region of small standard deviation (o> < 4) that extends 
across models within {R a ,i2iR^a) — (1-0,0.75) to (1.6,1.25) and defined with a slope (in 
rate multiplier ratios) close to unity with a spread of pa ±0.2 in R a ,i2- The IMF averaged 
production factors for the weak s-process isotopes are shown in Fig.|HJ The results for the 
individual masses can be found in the appendix for the weak s-only isotopes. 
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1.0 1.5 
,2 C(a,7) ,6 Rate Multiplier 



2.0 



Fig. 8. — Standard deviations for the IMF averaged production factors for the weak s-only 
isotopes, using the entire grid of initial masses. 

The results in Fig. [8] show a region of small standard deviation (o> < 4) that extends 
across models within (R a ,i2, R-3a) — (1-3, 1-0) to (1.5, 1.25) and defined with a slope close 
to unity with a spread of fh ±0.1 in R a ,i2- The production factors for all isotopes for the 
(Ra,i2, R3a) — (1.3,1.0) model are given in Fig.[9j This model lies within the region of 
minimum standard deviation for both the intermediate-mass and weak s-only isotopes. 
The neutrino-process isotopes, 7 Li, 11 B, 19 F, 138 La, and 180 Ta all show over-productions. 
Specifically, n B shows a produ ction factor very close to 16 {Pub/Pwo — 0.97), which 



agrees with 



Austin et al. 



(120111 ); these authors also show that this ratio varies by more 
than a factor of 2 at different values of R a ,i2- The low values for most of the heavy nuclei 
in Fig.|9]are expected, as we did not include the r-process or s-process contributions from 
AGB stars. 
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Fig. 9. — Production factors for all isotopes for the (R a ,u, R3a) = (1-3, 1.0) model. Shown 
are the 16 production factor (dashed line) with ±2 ranges (dotted lines). 



3.3. Implications for Stellar Remnants 

The analysis above assumed a perfect supernova success rate; i.e., abundances from all 
models were included even though some would collapse to a black hole without enriching 
the ISM with a SN event. We thus performed an additional analysis that removed the 
models that may result in possible "failed" supernovae, prior to IMF a veraging. Black 



hole f ormation following core collapse has been investigated recently by 



O'Connor and Qtt 



( 1201 lh . who identified a single parameter that can be used to roughly infer the fate of the 
core collapse event, using 62 progenitors. This compactness parameter is defined generally 
as, 
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( - M/Mg (4) 
?M R{M) /1000km' 1 J 

where M is the baryonic mass, and R (M) is the radial coordinate that encloses M at 
the time of core bounce. The releva nt specific £m for b l ack h ole formation is at a mass 
of M = 2.5 M , and £2.5 is used by lO'Connor and Ottl (1201 1[ ) to distinguish a possible 
boundary between successful and failed supernova explosions at the value of £ 2 .5 = 0.45. 
Their models with £2.5 < 0.45 were concluded to be likely successful supernova, using 
considerations of the time-averaged neutrino heating efficiency and subject (albeit mildly) 
to the equation of state (EOS) employed. 



A more recent analysis of over 100 supernova simulations by lUgliano et al.l (120121 ) have 
resultant NS and BH mass ranges that are compatible with a possible paucity of low mass 
BHs, which may imply a lower £ 95 value t han u sed by ?. A more refined boundary of 



£2.5 ~ 0.25 has been proposed by 



Woosleyl ( 120121 ). and will be adopted in the present work. 



In our analysis all models are assumed to have the same final kinetic energy for the 
ejecta (1.2 B), however, the explosion energies can vary. A larger explosion energy would 
cause successful SNe above £2.5 = 0.25, since now larger densities would be required to 
overcome this larger energy and prevent a successful SN explosion. It would also cause 
more material above the Fe core to escape the gravitational potential, resulting in a smaller 
mass cut and remnant mass for models already below this limit. 

The £2.5 values were computed for each model. The distribution of £ 2 .5 values for the 
25 M Q models is given as an example in Fig.HUJ Figures for the £2.5 values of the other 
masses can be found in the Appendix. 
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Fig. 10. — The distribution of £2.5 values for the 25 M models. 

For the 25 M Q models, those with low -R Qj i2 values are favored for successful SNe 
events, in addition to a local minimum in £ 2 .5 defined by a narrow strip close to the centroid 
value for this rate. In comparison, the 12 M , 13 M , 14 M , 15 M , and 16 M models 
all explode as successful SNe, whereas the 17 M models have only 6 of 176 that fail. The 
18 M models have 16 failed SNe, all above R a ,i2 > 0.9 and spanning the whole range of 
R^a. Among our models of 20 M , several fail towards high R Q) i2 and low R^ a values, with a 
fairly delineated boundary between successful and failed SNe beginning at a multiplier pair 
value of (^,12, #3a) = (1-2,0.75), and ending at (1.8,1.25). The 22 M and 27 M models 
show the same behavior as the 20 M models but with a boundary favoring more BHs, 
delineated by a slope defined by (1.0,0.75) and (1.4,1.25) for the former, and (0.5,0.75) 
and (0.8, 1.25) for the latter. Among our models of 30 M , there are 10 that explode at 
Ra,i2 = 0.5 along with a strip of successful SNe from (0.8, 0.75) to (1.3, 1.25). 
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The yields from the stellar winds of all models, as well as the explosive yields of the 
models that satisfy the condition £ 2 .5 < 0.25 were then averaged over an IMF. The result 
is given in Fig. [TT] for the intermediate-mass isotopes and in Fig. H2] for the weak s-only 
isotopes. 
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Fig. 11. — Standard deviations for the IMF- averaged production factors for the intermediate- 
mass isotopes. The yields from all stellar winds and the explosive yields from models that 
satisfy the condition £ 2 .5 < 0.25 were used in the averaging. 
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Fig. 12. — Standard deviations for the IMF-averaged production factors for the weak s-only 
isotopes. The yields from stellar winds and from the explosive yields of models that satisfy 
the condition £ 2 .5 < 0.25 were used in the averaging. 

As stated, Fig. [7| showed a small standard deviation region (op < 4). This region 
extended across models within {R a ,i2, R3a) — (1.0,0.75) to (1.5, 1.25) and was defined with 
a slope close to unity with a spread of ~ ±0.2 in R a ,i2- The impact of the compactness 
parameter for the intermediate-mass isotopes is that the (o> < 4) region from Fig.[7]is 
still observed in Fig. [HI but the latter has models in this region whose standard deviation 
has op larger by 1, although this effect may be less if an even finer mass grid was used. 
We conclude that the current 25 % uncertainty range of R a ,u agrees with observations. 
The region of small standard deviations for Fig. [IT] has a slope defined approximately by 
{R a , 12 ,R3a) = (1.1,0.75) to (1.6, 1.25) with a spread of « ±0.2 in R a<12 for R 3a < 0.95, and 
spread of Prf ±(0.3 — 0.4) in R a ,i2 for R 3a > 0.95. It is then possible to define the relation, 
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R aA2 = l-0R 3a + (0.25 ±0.3). (5) 

For a chosen i? 3cr value in the range 0.75 to 1.25, Equation^ gives a range of R a ^2 values 
that agree with observations, taking into account only the intermediate-mass isotopes. 

The impact of the compactness parameter on the weak s-only isotopes is that the 
(crp < 4) region from Fig. [8] becomes somewhat more sharply defined, in that the models 
above R 3a pa 0.95, from R a ,\2 ~ 1.0 to 1.3 have increased standard deviations by R3 1 — 2. 
In Fig. [T2l there is a region of small o~p with a slope defined by (1.1, 0.75) to (1.6, 1.25), with 
an average spread of ~ ±0.2 in R ay \2- This result should be interpreted with care, however, 
since the IMF averaging is subject to the w values, tak en from an approximat e analysis 



20121 ) which 



of weak s-process contributions to the solar abundances ( jWest and Hegerl 
uses only the relevant nuclear physics and does not employ stellar modeling. Furthermore, 
the decomposition of the weak s-process contributions does not address recent works 



that indicate possible evidence for an increase of t 



(IMashonkina et al. 



re s-elements in the Galactic disk 



2007 



Majorca et al. 



2011 



2012 



Jacobson and Friel 



20121 ). It is unclear 



what impact this would have on our analysis, however, since the weak s-isotopic abundances 
do not dominate their respective elemental abundances. Despite these issues, the weak 
s-process analysis in the present work is a step in the right direction; one simply cannot 
compare the weak s-process yields directly to the solar abundances (which contain main 
s-process abundances also). We thus caution the reader that whereas our weak s-process 
analysis is an improvement, it is also weakly constrained. 

The best values for the helium rates from our analysis should agree with observations 
for both the intermediate-mass and weak s-only isotope sets. We thus computed the average 
of the standard deviation values for both sets (Fig. [TBI. Whereas the optimal values for the 
helium burning rates from our analysis should reproduce the observed abundance ratios 
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for both the intermediate-mass and weak s-only isotope sets, the two sets will not be at 
the same level of production factor, because they have different astrophysical natures and 
galactic chemical evolution histories. For example, the weak s-process, being of secondary 
nature, is overproduced by a factor of m 2.0 at solar metallicity and less is made at lower 
metallicities, yet their relative abundance ratio should be about solar - and this is what we 
match. In contrast, the intermediate-mass isotopes are primary and should be produced 
at a solar level, and again within this subgroup the isotope ratios should be at solar 
level. Hence, they were analyzed separately before the results were combined, instead of 
computing the standard deviations for all isotopes as a single group. The ratio of weak 
s-process to intermediate-mass isotopes is not the solar abundance ratio and not expected 
to be, so fitting both at once would be wrong. 
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Fig. 13. — The average of the standard deviations for the production factors of the 
intermediate-mass (Fig. fTTj) and weak s-only isotopes (Fig. lT2"|) . 
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The region of small standard deviations for Fig. [13] has a slope defined approximately 
by (R aA2 ,R 3 a) = (1.1,0.75) to (1.6, 1.25) with a spread of « ±0.1 in R a>12 for R 3a < 0.95, 
and spread of ~ ±(0.2 — 0.3) in R aj ±2 for R 3a > 0.95. We then define the relation, 

R aA2 = 1.0R 3a + (0.35 ±0.2). (6) 

For a chosen R 3a value in the range 0.75 to 1.25, Equation[6] gives a range of R a ^i2 values 
that agree with observations, taking into account both the intermediate-mass and weak 
s-only isotopes. Note that this relation does not accurately define the regions of small 
standard deviation for either isotope set individually (see Fig.[TT]and Fig. n~2"j) . Note further 
that the region of small a> in Fig.[13]is "noisy," and has values along the line of best fit that 
are not minima of this region. We must re-emphasize that the analysis used is approximate, 
and the region contains values that may change if we change the parameters of the model 
or improve the completeness of the models included. Hence, whereas we can define this 
region, the analysis is likely insufficient to reliably distinguish between neighboring values 
within it. Subject to the approximations employed in the analysis, coordinate pairs for the 
reaction rate multipliers that satisfy this relation results in nucleosynthesis that equally 
agrees with current observations for both the intermediate and weak s-only isotopes, as far 
as the present study can determine. 

We also explored the dependence of our results on the IMF used in Equation[TJ 
We computed IMF averaged production factors for both the intermediate and weak s 
isotope list using two modified Salpeter IMFs (±0.3 and —0.3 added to the exponent). 
For the weak s-isotope list, both modified IMFs resulted in a difference in IMF-averaged 
production factors by < 6 % in the region of best fit rates identified by Equation^ For the 
intermediate-mass isotope list, both modified IMFs resulted in a difference in IMF averaged 
production factors by < 4 % in the region of best fit rates. We thus believe the choice of 
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IMF has only a small impact on the results, provided a reasonable IMF is chosen. 



The stellar model KEPLER, however, is only approximate: convection is treated 
using mixing length theory, effects of rotation, binary star evolution, and magnetic fields 
are ignored, many reaction rates and mass loss rates are not well enough known, the 
opacities have uncer t aintie s, etc. Some of these effects have been investigated. For example, 
Chieffi and Limongil (120121 ) studied the impact of rotation on solar metallicity massive stars 
in the range 13 — 120 M Q , and found over- productions o f F an d slight over-productions of 



weak s-process isotopes. Another study by llliadis et al. 
uncertainties that ir 
than those found by 



(120111 ) found several reaction rate 
uncertainties that in fluence massive star Al production, and found a range of 26 Al larger 



Tur et al. 



(120071 ). Although it was not the purpose of the present work 
to address the effects of the approximations in the stellar models, it is important to note 
that they can have a non-negligible impact on the results in some cases. Additionally, we 
interpolated yields across a finite IMF sampling and only for solar metallicity CCSNe stars 
instead of a full galactic chemical evolution model from big bang nucleosynthesis to the 
present Galaxy including all nucleosynthesis sources. For the isotopes we compare in this 
work, however, the assumption that solar composition stars should produce about their 
solar ratios is reasonable. 



3.4. Variations in Carbon Mass Fractions and Remnant Mass 



The baryonic mass of the progenitor of the remnant depends on the central carbon 
mass at the end of core-He burning^. An increase in R^ a or decrease in i? Q ,i2 results in an 
increase in the carbon abundance. An example of the resulting trend in baryonic mass is 



1 The gravitatio nal mass of the remn ant also depends on 



( iZhang et al.l 120081 1 . and equation of state ( iLattimer and Prakash 



;ype, formation scenario 



200 lh . 
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given in Fig. [HI for the 25 M Q models. Figures for the baryonic masses of the other models 
can be found in the Appendix. 
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Fig. 14. — The baryonic mass of the progenitor of the remnant for the 25 M Q models as a 
function of the R^ a and R a ,\2 multipliers. 

Fig. ITU shows that there is a decrease in baryonic remnant mass for increasing R^ a and 
decreasing R a ^2- This is because a larger carbon abundance at the end of core- He burning 
can support longer and more energetic carbon shell burnin g episodes, which allows the 



core to cool to lower entropy, yielding a smaller progenitor (jWoosley et al. 



2003 



Tur et al. 



20071 ). Note the apparent local maximum beginning at a R a ^2 multiplier value of ~ 0.7 
and extending to ~ 1.2, which coincides with the £2.5 local maximum in Fig.[10j This 
may be caused by non-convectiv e core-C burning t hat results in a more compact star and 



more massi ve baryonic remn ant (IHeger et al. 



observed by lTur et al. 



J2007j). 



20011 ). This non-monotonicity has also been 
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The cut-off mass for the remnant becoming a BH versus a NS is difficult to assess. A 
fraction of the NS progenitor is radiated away by neutrinos, which is dependent on the EOS 



(ILattimer and Prakash 



20011 ). A larger maximum for NS masses result from the "stiffest" 



EOS (those with larges t pressures 



has been calculated by 



Tolos et al. 



or a gi ven density), and an upper limit of Mhigh = 2.9 M Q 



(120121 ). Observational evidence places this limit between 



2.0 M Q < Mh ig h 2.5 M , although this may only be an indication of the limi t to the mass 



that can possibly be accreted in a binary system (ILattimer and Prakash 



20101 ) 



An example of the correlation between baryonic mass and central carbon mass at the 
end of core-He burning can be seen by comparing Fig.Owith Fig.lTSlfor the 25 M models. 
Figures for the carbon mass fractions of the other models can be found in the Appendix. 
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Fig. 15. — Central carbon mass fraction at the end of core-He burning for the 25 M models 
as a function of the R 3a and R a 12 multipliers. 



As shown in Fig.[15l high R 3a and low R a 12 multipliers result in higher carbon mass 
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fractions as expected, whereas low R 3a and high R ai u multipliers result in lower carbon 
mass fractions. Comparing Fig. and Fig. IT51 shows an overall inverse relation between the 
carbon mass fraction and the remnant mass (see also Fig.[TH]and Fig. l2"0"j) . 



4. Conclusions 

This paper studies the effect of changing the helium burning rates 12 C(a,7) 16 and 
4 He(2a, 7) 12 C independently to map the effects on stellar evolution and nucleosynthesis. We 
follow the entire evolution from hydrogen burning through the SN explosion, including wind 
and SN yields, and considering fallback, mixing, and which stars make SNe (£2.5 < 0.25) or 
collapse to BHs without SN, and finally integrating the yields over an IMF. In total, we 
calculated nucleosynthesis for a grid of 176 models for each of 12 stellar masses from 12 M 
to 30 M (2112 models). This is by far the most extensive investigation on the effects of 
rate variatio ns for the heliu m burning reactions to date, and the first to use updated solar 



abundances ( jLodders 



20091 ). 



Combining constraints on intermediate-mass and weak s-only isotopes, we find a best 
fit for rate multipliers (R a ,u,R3a) = (1.1,0.75) to (1.6,1.25) with an average spread of 
« ±0.2 in R a ,i2- More generally, we find a relation between R a ,i2 and Rs a for good fits to 
nucleosynthesis given by R a: n = + (0.35 ± 0.2) in the range R 3a = 0.75 to 1.25. 

We also provide the line of best fit using only the intermediate- mass isotope list, given by 
^,12 = 1.0^ + (0.25 ±0.3). 

In this analysis, all models are assumed to have a final kinetic energy of the ejecta 
of 1.2 B. Real supernovae have a range of explosion energies. A larger explosion energy 
would allow successful SNe above the £2.5 = 0.25 limit, since now larger densities would be 
required to overcome this larger energy to prevent a successful SN explosion. It would also 
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cause more material above the Fe core to escape the gravitational potential, resulting in a 
smaller mass cut and remnant mass for models already below £ 2 .5 = 0.25. 

Small changes in the reaction rates can result in significant differences in the convection 
structure of the star, which is not just a numerical artifact but due to the physics of shell 
burning. This introduces "noise" into the comparisons, and is a warning for calculations 
done for specific cases; the general conclusion may be influenced by isolated processes. 
Given an astrophysical model that includes all the important physics with perfectly known 
physical input parameters, the values for the standard deviations shown in the figures 
should have a minimum region reflecting only the uncertainties in isotopic abundances. 
We know, however, that the stellar model used is only approximate and can affect the 
nucleosynthesis, as discussed in Section 3.3. For the isotopes we compare in this work, 
however, the assumption that solar composition stars should produce about their solar 
ratios is reasonable. 

If we change the parameters of the models or improve the completeness of the models, 
we expect the best fit rates to change also. Hence it is not necessarily the case that the best 
fit rate for an observable (in our case the abundances) coincides with the true rate. This 
suggests that the best reaction rates we obtain in our analyses are at some level effective 
rates. Analogous procedures have been used in many areas of physics. For example, in 
the shell model of nuclear physics, an effective nucleon-nucleon interaction (close to but 
not identical to the true interaction) is chosen to fit the low-lying spectra of many nuclei, 
and this reaction is very successfully used to predict other observables. Similarly, we think 
that the present procedure can provide better predictions of other astrophysical quantities, 
remnant masses or neutrinos synthesis of certain isotopes, for example. 

Whereas our comparisons do not have the power to accurately determine the helium 
burning reaction rates, they do show that the experimental values are not too far from the 
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truth and that changes needed to compensate (in an effective interaction sense) for model 
uncertainties are not large. Thus, similar calculations necessary to assess the situation, 
as overall model uncertainties decrease, will be less demanding; we have shown that a 
significant part of the uncertainty space is irrelevant. It is also clear that if one of the two 
helium burning reactions is much better determined, the effective rate for the other will be 
much better determined. For example, if it is later determined that is near 1.25 times 
the present experimental centroid value, then the best fit R a ,i2 will be 33 % larger than 
the experimental centroid value. It may also happen that if both reactions become well 
determined, they would not agree with the effective interaction that best reproduces the 
abundances. In such an event, this work may still serve to provide an evaluation of other 
model uncertainties and point the way to improvements. 
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A. Appendix 
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Fig. 16. — Standard deviations of the production factors of intermediate-mass isotopes for 
all models as a function of the R^ a and R a ,u multipliers. The horizontal axis gives the R a ,i2 
multiplier, the vertical axis gives the Rz a multiplier, and the color scale gives the logarithm 
of the standard deviations. 
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Fig. 17. — Standard deviations of the production factors of weak s-only isotopes for all 
models as a function of the i? 3Q and -R Qi i2 multipliers. The figure follows the convention of 
Fig.dS 
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Fig. 18. — Compactness parameter values for all models as a function of the R?, a and R a ^i2 
multipliers. The color scale gives the logarithm of the compactness parameter, and the rest 
of the figure follows the convention of Fig.lT6l 
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Fig. 19. — Baryonic mass of the progenitors (M ) for all models as a function of the R^ a 
and -R Qj i2 multipliers. The color scale gives the logarithm of the baryonic mass, and the rest 
of the figure follows the convention of Fig.lT6l 
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Fig. 20. — Central carbon mass fractions at the end of core-He burning for all models as a 
function of the R 3a and -R a ,i2 multipliers. The color scale gives the carbon mass fraction, 
and the rest of the figure follows the convention of Fig.[T6l 
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